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Abstract 

We discuss the modification of a jet fragmentation function due to medium-induced 
partonic energy loss in context of leading particle observables in ultrarelativistic nucleus- 
nucleus interactions. We also analyze the relation between in-medium softening jet frag- 
mentation function and suppression of the jet rates due to energy loss outside the jet cone. 
The predicted anti-correlation between two effects allows to probe a fraction of partonic 
energy loss carried out of the jet cone and truly lost to the jet. 
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1 Introduction 



Jets play one of the central roles as a promising tool to study properties of quark-gluon plasma 
(QGP) expected to be created in heavy ion collisions at RHIC and LHC. Medium-induced 
energy loss of energetic partons, the so-called jet quenching, has been proposed to be very 
different in cold nuclear matter and in QGP, resulting in many challenging observable phenom- 
ena (see for the review [1] and references therein). In particular, softening jet fragmentation 
function and, as a consequence, suppression of high-pr hadron spectrum in nucleus-nucleus 
collisions relative to their production in independent nucleon-nucleon interactions, are consid- 
ered [2, 3, 4, 5, 6, 7, 8]. Recent RHIC data on inclusive high-p^ charge and neutral hadron 
production from STAR [9], PHENIX [10], PHOBOS [11] and BRAHMS [12] experiments show 
such kind of suppression and are in agreement with the jet quenching hypothesis. However, 
since at the moment direct event-by-event reconstruction of jets and their characteristics is not 
available in RHIC experiments, the assumption that integrated yield of all high-p^^ particles 
originates only from jet fragmentation is not fully clear (see, for example, [13]). At the LHC, a 
new regime is reached where hard and semi-hard QCD multi-particle production can certainly 
dominate over underlying soft events [14]. CMS experiment at LHC [15] will be able to provide 
adequate jet reconstruction using calorimetric measurements [16, 14]. Thus identification of 
leading particle in a jet (i.e. particle carrying the maximal fraction of jet transverse momen- 
tum) allows the measurement of jet fragmentation function (JFF) to be done. Comparison of 
JFF in AA and pp collisions (or in central and peripheral AA interactions) may give information 
about in-medium modification of JFF. 

The crucial related question here is: how much energy loss falls outside the typical jet cone 
and is truly lost to the jet? There are some discussions on this subject in the literature [17, 
18, 19, 20, 21]. In fact, since coherent Landau- Pomeranchuk-Migdal radiation induces a strong 
dependence of the radiative energy loss of a jet on the angular cone size, it will soften particle 
energy distributions inside the jet, increase the multiplicity of secondary particles, and to a 
lesser degree, affect the total jet energy. On the other hand, collisional energy loss turns out 
to be practically independent of jet cone size and causes the loss of total jet energy, because 
the bulk of "thermal" particles knocked out of the dense matter by elastic scatterings fly away 
in almost transverse direction relative to the jet axis [17]. Moreover, the total energy loss of a 
jet will be sensitive to the experimental capabilities for low-p^ particles, products of soft gluon 
fragmentation. For example, in the CMS case, most of these low-py hadrons may be cleared 
out of the central calorimeters by the strong magnetic field [15, 16]. 
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In this Letter we analyze relation between in-medium softening JFF and suppression of jet 
spectrum due to energy loss outside the jet cone. 



2 Medium-modified jet fragmentation 

Let us recall that in the leading order of perturbative QCD the jet production cross section 
with transverse momentum pt and rapidity y in a nucleon-nucleon collision is given by 

where Xi — Pi/pa, xj — Pj/pb are the initial momentum fractions of nucleons a, b carried by the 
interacting partons of types i, j; fl{xi, Q^), fl{xj, Q^) are the parton distribution functions for 
the colliding nucleons a, b; s, t and u are the Mandelstam variables of hard parton subprocess; 
da/dt{ij — > kl) is the Born cross section for the hard ij kl scattering subprocess; = 

p"^ = {tu)/s; t/u = (xi/xj) exp (— 2y); ^ = pr/ \/s{Gxp {y)/xi + exp {—y)/xj) is the momentum 



fraction of parton k carried by the jet; y/s = ys/{xiXj) is the center of mass energy of the 
colliding nucleons. To calculate the inclusive cross section for "jet-induced" hadron production 
with transverse momentum pj^ and rapidity y'^ in a nucleon-nucleon collision, one convolutes 
the jet production cross section with the fragmentation function D^{z',p^) for the parton of 
type k into hadron h: 

-m^'yr), (2) 



d{p!},)'^dy^dz' dp^dy z'"^' 
where y^ — y, and z' — p'^/p'' — p'^/px is the momentum fraction of a parton k carried by the 
final observable hadron. 

The parton distribution functions fl{xi, Q^) are measured in deep inelastic scattering ex- 
periments such as those at HERA [22], while fragmentation functions D^{z', Q"^) are extracted 
from e^e" annihilation from PETRA, PEP and LEP [23] and from hadronic collisions from 
UAl [24] (gluon fragmentation function). 

In nuclear interactions medium-induced energy loss of fast partons can modify the cross 
sections for high-p^ hadrons and jets together with other potentially important nuclear effects 
like, for instance, parton shadowing playing significant role at small values of the momentum 
fraction carried by the interacting partons. However for sufficiently hard jets and hadrons 
under consideration [xi^j > 0.2) this effect is neghgible as well as next-to-leading order (NLO) 
corrections {K factor ~ 1) [25]. As a result the rate of k-type jets of finite angular cone size 
6o in mid-rapidity with transverse momentum pj^ in AA collisions at given impact parameter 
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b can be estimated as 



db^^^"^' ^ / / ^^^^^(^^^^-^^^^^ ^^P^ ' 

where ri^2ib,r,ip) are the distances between the nucleus centers and the jet production vertex 
y(r cos'^, rsin'^); rmaxi^ii^) < Ra is the maximum possible transverse distance r from the 
nuclear coUision axis to the V; Ra is the radius of the nucleus A; TA{ri^2) is the nuclear 
thickness function (see Ref. [25] for detailed nuclear geometry explanations). The effective shift 
Ap^*(r, "0, ^o) of jet momentum spectrum depends on the jet cone size ^o- 

The partons will not hadronize inside QGP. In a hadronic medium, we assume that the 
fragmentation functions can be approximated by their forms in vacuum (see, however, Refs. [26, 
27] where the corrections to the parton fragmentation function in a thermal medium were 
discussed) . We take into consideration the fragmentation of leading partons only and omit the 
fragmentation of emitted gluons, because we are interested here in the leading particle in a jet. 
Then the rate of high-p^- jet-induced hadrons can be estimated as 

d{pl^)Hydz'^^^ ^ I "^^^^^^^^^^(^^^ d^p^Ty Y^^^^' '^^^ ' 

where the shift Ap^ of hadron momentum distribution generally is not equal to the mean 
in-medium partonic energy loss due to the steep fall-off of the p^-spectrum [3] . 

The integral jet suppression factor Q can be introduced by the natural way as the ratio of 
jet rate with energy loss to jet rate without one, 

Q-W'„,J^ / / .«Y*;^(A# ^ 0). (5) 



d{t?T ydyl J d{j?T fdy 

Pt min Pt min 



JFF is defined here as 



/ ,m^ay,,^J^,s(.-/,ir>?)/ j (6) 

^'Pt min I^T min 

and coincides approximately with D^{z, (Pr*min)^) ^^r the case without energy loss and with the 
type of jet specified. Note that z = Pt/Pt^{— z'Pt/Pt^) is experimentally observable quantity 
depending on jet cone size ^o- Further we are interested also in the ratio 

1 



(7) 
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where the integral hadron suppression factor 



ZO Jot 

* Pt min 



^ fr min 



d{p^ydydz' 
d{pj^ydydz' 



r c dN^^^^ 

J,. I rf(pJ)V,rf/ -_(Ap, = 0)^(.-p^/p¥') (8) 



is distinguished from the differential hadron quenching factor usually defined as [3, 6, 7, 8] 



3 Model 



In order to generate the initial jet distributions in nucleon-nucleon sub-collisions at ^/s = 
5.5 TeV, we have used PYTH1A_5.7 [28]. After that we perform event-by-event Monte-Carlo 
simulation of rescattering and energy loss of partons in QGP (for model details one can refer 
to our previous papers [17, 25]). The approach relies on an accumulative energy losses, when 
gluon radiation is associated with each scattering in expanding medium together including the 
interference effect by the modified radiation spectrum as a function of decreasing temperature 
dE/dl{T). The basic kinetic integral equation for the energy loss AE as a function of initial 
energy E and path length L has the form 

^ ^dP{l) ^^^^dE{l,E) dP{l) _J_ 

A(I) 

where I is the current transverse coordinate of a parton, dP/dl is the scattering probability 
density, dE/dl is the energy loss per unit length, A = Xjiap) is in-medium mean free path, 
p oc is medium density at temperature T, a is the integral cross section of parton interaction 
in the medium. Such numerical simulation of free path of a hard jet in QGP allows us to obtain 

any kinematical characteristic distributions of jets in final state. Besides the different scenarios 
of medium evolution can be considered. 

For the calculations we have used coUisional part of loss [25] , 



3TE/2 



dl ATXa 



and the dominant contribution to the differential cross section 

da ^ 27ral(t) Un 

dt ' ' {?,?,- 2Nf)\n{t/KlcD) ^ ^ 
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for scattering of a parton with energy E off the "thermal" partons with energy (or effective 
mass) mo ~ ST ^ E. Here C = 9/4,1,4/9 for gg, gq and qq scatterings respectively, 
is the QCD running coupling constant for Nf active quark flavours, and Aqcd is the QCD 
scale parameter which is of the order of the critical temperature, A.qcd — Tc 200 MeV. The 
integrated cross section a is regularized by the Debye screening mass squared fij^iT). 

The energy spectrum of coherent medium-induced gluon radiation and the corresponding 
dominated part of radiative energy loss was estimated using BDMS formalism [18] : 



where ti = L/ (2Ag), y = uj/E is the fraction of the hard parton energy carried by the radiated 
gluon, and Cr — 4/3 is the quark colour factor. A similar expression for the gluon jet can be 
obtained by substituting Cr — 3 and a proper change of the factor in the square bracket in (13), 
see Ref. [18]. The integral (13) is carried out over all energies from (^min = Elpm — {\ 
is the gluon mean free path), the minimal radiated gluon energy in the coherent LPM regime, 
up to initial jet energy E. Note that although the radiative energy loss of an energetic parton 
dominates over the coUisional loss by up to an order of magnitude, the relative contribution of 
collisional loss of a jet growths with increasing jet cone size due to essentially different angular 
structure of loss for two mechanisms [17]. 

The medium was treated as a boost-invariant longitudinally expanding quark-gluon fluid, 
and partons as being produced on a hyper-surface of equal proper times r [29]. For certainty 
we used the initial conditions for the gluon-dominated plasma formation expected for central 
Pb-Pb coUisions at LHC [30]: Tq ~ 0.1 fm/c, Tq ~ 1 GeV, Pg ^ 1.95T^. For non-central 
collisions we suggest the proportionality of the initial energy density to the ratio of nuclear 
overlap function and effective transverse area of nuclear overlapping [25]. 

Thus we consider that in each i-th scattering off the comoving particle (with the same lon- 
gitudinal rapidity) a fast parton loses energy collisionally and radiatively, Ae^ = tj/(2mo) +uji, 
where and uji are simulated according to Eqs. (11) and (13) respectively, and the distribution 
of jet production vertex - according to Eq. (3). Finally we suppose that in every event the en- 
ergy of an initial parton decreases by value ApT{r, ip) = Y.i and the jet energy loss, Apj^, is 
the product of e times energy loss of an initial parton, Ap^^*(r, i/j^Oo) — e- Apr(r, -0). We exam- 
ine here the fraction e of partonic energy loss carried out of the jet cone as a phenomenological 
parameter, since the treatment of angular spectrum of emitted gluons is rather sophisticated 




(13) 



(14) 
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and model-dependent [17, 18, 19, 20, 21] ^. We was interested in relatively realistic values of e 
in the range from to 1 (although e can be even larger than 1 at small ^o, see [18, 20]). 

In our simulations we have considered channel with neutral pions only as leading particles 
in jets. One can expect the similar results for leading charged hadrons. Prom methodical point 
of view tracking in heavy ion environment at LHC is rather complicated (although solvable) 
task, while the reconstruction of electromagnetic clusters in calorimeters being at the moment 
more understandable [14, 16]. At high enough transverse momentum of ttq (> 15 GeV at CMS 
case [14]), two photons from pion decay fall into one crystal of electromagnetic calorimeter, and 
traditional technique for reconstructing ttq's using two-photon invariant mass spectrum does 
not work. However, such electromagnetic cluster can be identified as a leading ttq, if it belongs 
a hard jet and carries the significant part of jet transverse energy. 

4 Numerical results and conclusions 

Figure 1 shows JFF (6) (without the jet-type specification and therefore experimentally observ- 
able) for leading tto's for the cases without and with medium-induced energy loss obtained in the 
frameworks of our model [17, 25] in central and minimum-bias Pb— Pb collisions. The threshold 
for jet reconstruction P7^\nin = 100 GeV was used [16] ^. If £ close to ( "small- angular radia- 
tion" dominates), then the factor of jet suppression Q^^^ (5) is close to 1 (there is almost no jet 
rate suppression), and effect on JFF softening is maximal. Increasing e (the contribution from 
"wide-angular radiation" and collisional loss grows) results in stronger jet suppression {Q^^^ 
value decreases), but effect on JFF softening becomes smaller, especially for highest z. Figure 
2 presents the £-dependences for jet suppression factor Q^^* (without the jet-type specification) 
and ratio (7) of JFF with energy loss to JFF without loss, D{z > zo)/D{z > zq, Apr — 0), for 
zq — 0.5 and 0.7. Note that in the case without jet quenching, fraction of events when leading 
TTo carries larger 50% (70%) of jet transverse momentum is 6.3 • 10~^ (9.3 • 10~^). One can see 
the distinctive anti-correlation between strengthening jet suppression and JFF softening (ratio 
(7) can be even greater than 1 at large enough e and z values) determined mainly by fraction 
of partonic energy loss carried outside the jet cone. The physical reason for the effect to be 

^Of course, this parameter is regulated by the jet cone size 9o and its dependence on 9o was numerically 
investigated [18, 20]. We can formally recalculate our result as a function of using such numerical estimation 
or taking into consideration the angular structure of radiative and collisional energy loss in each scattering com- 
plicating our simulation. However this analysis demands also the jet finding procedure specification and other 
technical details. We believe that our simplified treatment here is enough to demonstrate the anti-correlation 
between jet suppression and JFF softening, moreover the integral hadron suppression factor <3''(PT*inin' -^o) in 
Eq. (7) depends on s implicitly (only via z-definition). 

^The estimated statistics for > 100 GeV jets in CMS acceptance is high enough, at the level ~ 10^ jet pairs 
per 1 month of LHC Pb-Pb run [16, 14] 
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opposite in the jet suppression factor and the fragmentation function is it follows. Increasing 
e results in decreasing final jet transverse momentum, p!^^ — p^^*(ApJ^* — 0) — e • Ap^^* (which 
is the denominator in definition of 2; = p^/Pr* '^^^ (6)) without any influence on the nu- 
merator of z and, as a consequence, in reducing effect on JFF softening, while the integral jet 
suppression factor (5) becomes larger. The remarkable prediction here is that the effect on 
jet rate suppression becomes comparable with the effect on JFF softening at quite reasonable 
value s ~ 0.3. 

In summary, we have analyzed the relation between in-medium softening jet fragmentation 
function (in leading ttq channel) and suppression of the jet spectrum due to energy loss outside 
the jet cone. We beheve that this kind of analysis can be performed for the heavy ion coUisions 
at LHC experiments. The observation of significant JFF softening without substantial jet rate 
suppression would be an indication of the fact that small-angular gluon radiation is dominating 
mechanism for medium-induced partonic energy loss. Increasing contribution of wide-angular 
gluon radiation and coUisional loss can result in jet rate suppression, while the effect on JFF 
softening becomes less in this case. If the contribution of the "out-of-cone" jet energy loss is 
large enough, the jet rate suppression may be even more significant than JFF softening. 

Acknowledgements. Discussions with Yu.L. Dokshitzer, L.I. Sarycheva, I.N. Vardanian, 
U. Wiedemann, P. Yepes, B. Zakharov and G.M. Zinovjev are gratefully acknowledged. 
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z=PT(^o)/PT(jet) 



Figure 1: Jet fragmentation function for leading ttq's without (solid curves) and with medium- 
induced energy loss for £ = (dashed curves), 0.3 (dash-dotted curves) and 0.7 (dotted curves) 
in central (a) and minimum- bias (6) Pb-Pb collisions. > 100 GeV and | |/^"* |< 3. 
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Figure 2: Jet suppression factor Q^'^^ (solid curves) and ratio of JFF with energy loss to JFF 
without loss, D{z > zq)/ D{z > zq, Apr = 0), for zq = 0.5 (dashed curves) and 0.7 (dash-dotted 
curves) in central (a) and minimum- bias (6) Pb-Pb collisions. > 100 GeV and | |/^"* |< 3. 
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